Transport Properties of a spinon Fermi surface coupled to a U(l) gauge field 



o 
o 

W)' 
< 



c3 



i 

c 

o 
o 



> 
o 
in 
oo 

00 

o 

O 



Cody P. Nave and Patrick A. 
Department of Physics, Massachusetts Institute of Technology, 

(Dated: February 1, 2008) 



Lee 

Cambridge, Massachusetts 02139 



With the organic compound k-(BEDT-TTF)2-Cu2(CN)3 in mind, we consider a spin liquid system 
where a spinon Fermi surface is coupled to a U(l) gauge field. Using the non-equilibrium Green's 
function formalism, we derive the Quantum Boltzmann Equation (QBE) for this system. In this 
system, however, one cannot a priori assume the existence of Landau quasiparticles. We show that 
even without this assumption one can still derive a linearized equation for a generalized distribution 
function. We show that the divergence of the effective mass and of the finite temperature self- 
energy do not enter these transport coefficients and thus they are well-defined. Moreover, using a 
variational method, we calculate the temperature dependence of the spin resistivity and thermal 
conductivity of this system. 



PACS numbers: 74.70.Kn, 71.10.Hf, 71.18.+y 



I. INTRODUCTION 



Recent experiments have shown evidence that the or- 
ganic compound k-(BEDT-TTF) 2 -Cu 2 (CN) 3 maybe the 
first experimentally realized spin liquid in dimension 
greater than onei This quasi two-dimensional material 
can be described as an effectively isotropic spin 1/2 sys- 
tem on a triangular lattice. It is found experimentally to 
be insulating and shows no evidence of long range mag- 
netic order down to mK temperatures. Fitting the sus- 
ceptibility using the high temperature series expansion 
of the spin 1/2 Heisenberg model on a triangular lattice, 
the exchange coupling J is roughly 250A". The static 
spin susceptibility also remains finite down to the lowest 
temperatures measured. 1 Because of the lack of magnetic 
order even at temperatures many orders of magnitude 
lower than the exchange coupling J and the experimen- 
tal evidence for abundant low energy spin excitations, 
there has been a proposal that this system may be well 
described by a spin-liquid where a spinon Fermi surface 
is coupled to a U(l) gauge fieldi&si Recent work with this 
model has lead to the suggestion of a possible spinon 
pairing state that may explain the observed features in 
the experimental measurements of the specific heat and 
magnetic susceptibility^^ In this paper, we focus on the 
original model, i.e. we consider temperatures above the 
pairing transition temperature but small compared to the 
exchange temperature. In this regime, we start with our 
model Lagrangian and proceed to develop a version of the 
quantum Boltzmann equation (QBE). Despite potential 
pitfalls that we discuss below, we show that the QBE 
is well-defined and apply it to systems with steady-state 
thermal and spin currents. We show that the transport 
coefficients are finite and in particular calculate the tem- 
perature dependence of the spin resistivity and the more 
experimentally accessible thermal conductivity. 

We begin by considering the t-J model on the trian- 
gular lattice. We construct a mean-field state by ap- 
plying the slave-boson formalism and enforcing the local 
constraint of no double occupancy exactly. It is known 
that considering fluctuations around this mean-field state 
leads to a U(l) gauge theory^ More recently, the slave- 
rotor representation has been applied to the Hubbard 



model on a triangular lattice and it was shown that again 
one can arrive at a U(l) gauge theory.^ Because of the 
large number of low energy excitations due to presence 
of the spinon Fermi surface, we assume that a deconfincd 
state occurs and thus we consider a non-compact U(l) 
gauge theory. In other words, we assume that instanton 
effects are negligible. Thus our starting point is the La- 
grangian for a 2-D spinon Fermi surface system coupled 
to a non-compact U(l) gauge field given by 

1 2 

C = tjj* (d - ia - fjt) Vv + tt— C (~ iV ~ a ) Vv, (1) 

Zm 

where the gauge field kinetic energy term has been 
dropped because its strength is inversely proportional to 
the charge gap which is large since we are assuming we 
are in the insulating phase. tp a is the spinon field and 
the gauge field is a = (ao, a), fi is the chemical potential. 
We work in Coulomb gauge V • a = 0. 

We can then proceed to integrate out the spinons in or- 
der to generate dynamics for the gauge field. This cannot 
be done exactly; however, we can work in the Gaussian 
approximation also known as the random-phase approx- 
imation (RPA) . We then consider spinons coupled to the 
effective action for the gauge field coming from the spinon 
bubbles. The use of the RPA can be formally arranged 
in the standard way through the 1/JV expansion by in- 
troducing N species of fermions. 8 

The longitudinal part of the gauge propagator is re- 
lated to the density-density response and thus does not 
show any singular behavior for low energies and momen- 
tum. The transverse part however does give rise to long 
range interactions. After integration, we find that the ef- 
fective action for the gauge field is S(a) — J2 q U(q)a q a q , 
where 
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V^q 2 + q$ + q 
and 7 = — . Thus the effective gauge 



(2) 



where xb 



field propagator is given by D(q) — n(g) _1 . Rotating 
back to real time, qo — iv, and working in the limit that 
v -c fF<7> the gauge propagator becomes 



D(q,v) = 



1 



2 " 



(3) 



2 



where q is now the magnitude of q£>2 

For the remainder of the paper we consider the effec- 
tive theory given by taking the Lagrangian of Eq. [1] and 
adding gauge field dynamics through the RPA propa- 
gator of Eq. [3l This particular gauge theory has been 
studied previously in the context of the half-filled Landau 
level (Ref. IX Op and high-temperature superconductors 
(Refs. and llll) In particular, the spinon self-energy 
correction due to the RPA gauge propagator has been 
examined. Again one finds that the most singular cor- 
rection comes from considering the transverse gauge field 
fluctuations. To one-loop order, the self-energy S(k, w) 
is found to be ReS ~ ImS ~ w 2 / 3 . We note that this 
implies a vanishing quasiparticle spectral weight, i.e. the 
Landau criterion for quasiparticles is violated and thus 
the Fermi liquid picture is invalid for this system. More- 
over, the effective mass is found to diverge at the Fermi 
surface as £^ where £k = £ k — M = k 2 /2m — /z.— 

Despite the fact that quasiparticles are ill-defined in 
this system, we examine the standard expressions for the 
spin resistivity and thermal conductivity. The spin re- 
sistivity is given by ps = — where 1/r is the momen- 
tum relaxation rate. From the self-energy correction to 
the fermion propagator, we calculate the momentum re- 
laxation rate 1/r ~ T~ 4 / 3 £ Beyond the assumption of 
the validity of the quasiparticle picture, in order to ar- 
rive at ps, we also need to consider the effective rather 
than bare mass for the spinons; however, as mentioned 
above the effective mass is divergent. In section IV, we 
find that p$ ~ T 4 / 3 , a result which is consistent with 
the calculation when the quasiparticles are assumed to 
be well-defined and the renormalization of the mass is 
ignored. 

The standard simple result for the thermal conduc- 
tivity gives k ~ Cv 2 te where G is the specific heat, 
v is the particle velocity and 1/te is the energy relax- 
ation rate. Again from the self-energy correction to the 
spinon propagator, we calculate the energy relaxation 
rate 1 /rg ~ T 2 / 3 £ For a system of fermions C — 7T and 
the velocity would be temperature independent. Thus 
k/T ~ T~ 2 / 3 . However these assumptions are again not 
justified in view of the divergent effective mass and lack of 
well-defined quasiparticles. Considering the mass renor- 
malization leads to a specific heat G ~ T 2 / 3 and thus 
k/T ~ T -1 ; moreover, the velocity goes to zero as the ef- 
fective mass diverges. Thus it is unclear how to proceed. 
In section V, we see that the power law T dependence of 
the thermal conductivity given by the naive arguments 
ignoring the effects of the effective mass turns out to be 
correct. 

Because of these issues, we are forced to consider the 
interactions between the spinons and gauge bosons more 
carefully and thus turn to a quantum Boltzmann de- 
scription of the system. As mentioned above, we can- 
not derive a QBE using the Landau quasiparticle picture 
because Fermi-liquid theory is invalid for this system. 
We find however that we can proceed by following the 
work of Prange and Kadanoff who studied the electron- 
phonon system at temperatures high compared to the De- 
bye temperature* 1 ^ At high temperatures, the electrons 



rapidly emit phonons so that their precise energy is not 
well defined. Thus they deal with an analogous situation 
where the quasiparticle picture breaks down. Closely fol- 
lowing their work, we find that if the self-energy at small 
frequencies is independent of £k, we can define a gener- 
alized distribution function and derive a closed equation 
describing the dynamics of this generalized distribution 
function. This equation is analogous to the standard 
Fermi- liquid QBE. We note that the derivation of the 
QBE for this system is very closely related to the work 
done by Kim et al. in Ref. [H studying the v = 1/2 
fraction quantum Hall state, except that we derive the 
QBE in a different coordinate system and linearize in a 
different way, which allows us to use variational methods 
to calculate the transport properties. Note that a sim- 
ilar derivation of the QBE for generalized distribution 
functions is also considered in the work by Mahan.— 



II. DERIVING THE QBE 

To derive the QBE for this system, we work in the 
standard non-equilibrium Green's function formulation. 
We begin with two matrices of Green's functions G and 
£ that satisfy Dyson's equation 



G — GqYjG, 



where 



G = 



G> 



-G< 



(4) 



(5) 



with £ defined similarly. Note that following Ref. 111 
Eq. [¥] the product of two functions actually implies an 
integration over a shared space-time variable. Also here 
we use the standard definitions, following Refs. [Ill and 

EE 



G > (xi,x 2 ) 
G < {xi,x 2 ) 
G t {xi,x 2 ) 

Gi{xi,x 2 ) 



-i (ip{xi)^ (x 2 )) 
i (x 2 )ip(xi)) 

e(t 1 -t 2 )G > {x 1 ,x 2 )^ 
e(t 2 -t 1 )G < (x 1 ,x 2 ) 
e(i 2 -ii)G > (x 1 ,x 2 )- 
e(ti-i 2 )G < (x 1 ,x 2 ), 



(6) 
(7) 

(8) 
(9) 



with associated self energies S > , S t and £j. Here 
x = (r,t). Go denotes the non-interacting Green's func- 
tions. These Green's functions are related to the stan- 
dard retarded (G R ) and advanced (G A ) Green's functions 
through 



G 



Gt — G' 



G< - Gt 



(10) 
(11) 



We perform a change of variables so that all the Green's 
functions are expressed in terms of relative and center of 
mass like coordinates. Throughout this paper, we work 
with the Fourier transform of the relative coordinates so 
that we can write G < (k, ui, r, t). 
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For a general system of fermions in thermal equilib- 
rium, we can write that 



G<(k,w) = if (uj)A(k,oj) 
G>(k,uj) = -i(l-/ (w))yl(k,w) 



(12) 
(13) 



where fo(u)) is the Fermi distribution function at some 
temperature T. A(k,w) = -i{G R {k,uj) - G A (k,uj)) is 
the spectral function given by 



A(k,u) = 



-2lmY, R (k,uj) 



Re£«(k,w)] 2 + (ImE«(k,w)) 2 ' 

(14) 



Kadanoff used in Ref. [l2j to derive a generalized QBE for 
the electron-phonon system that is valid at temperatures 
high relative to the Debye temperature. 

Following the work of Ref. [12, we assume that the 
system remains close enough to local equilibrium that the 
self-energy is independent of £ at all times. Combining 
this assumption with the fact that the J -^-A — 1, it 
follows that G < and G > are sharply peaked functions of 
£. Integrating over the region of peaking in £, we can then 
define the generalized distribution function /(k, to, r, t) as 



2tt 



-iG<^,k,uj,T,t) =/(k,w,r,t), 



(16) 



In Fermi liquid theory, the quasiparticles are well de- 



fined because Im H R 



<C uj for small uj. This means 



that the equilibrium spectral function is sharply peaked 
as a function of ui, so that ignoring the incoherent back- 
ground, it can be written as 

A(k, uj) = 2ir5 (u-Zk- Re£ fl (k, uj)) . (15) 

In this paper we are considering linear response, so we 
assume that the system is slowly varying in space and 
time and that therefore there is a notion of a local equi- 
librium temperature T for every (r, t). Assuming that 
the system remains close enough to equilibrium that the 
sharp to peaking of the spectral weight remains valid, the 
standard Landau quasiparticlc QBE for the fcrmion dis- 
tribution function f(k, r, t) then follows. 

As mentioned in the Introduction, in this model which 
is described by Eqs. [T]and[31 both the real and the imag- 
inary parts of the self-energy of the spinons scale as lu 2 ^ 3 
for small uj. This violation of the Landau criterion for 
the existence of well-defined quasiparticles invalidates the 
normal derivation of the QBE since the spectral weight is 
no longer sharply peaked in uj and thus cannot be written 
in the form of Eq. [151 

We find that we can still proceed to derive a QBE for 
this system due to the form of the self-energy. First we 
change variables from k to £ = £k and k. Then we note 
that since the self-energy is independent of the magnitude 
of k, we can write E fl (k, uj) — £"(k, cj), and that then for 
small enough uj, A(k, uj) is a peaked function of £ around 
£ = 0. This property of the invariance of the self-energy 
with the magnitude of k, i.e. that it is only a function of 
uj and k, is exactly the same property that Prange and 



J 



where /(k, lu, r, t) is the density of spinons with momen- 
tum in the k direction, energy uj at a given position r 
and time t. Similarly we have 



rf£ 
2tt 



*G>K,£,w,r,t) =l-/(k,w,r,t). (17) 



Having established the definition of the generalized dis- 
tribution function in a system without well-defined Lan- 
dau quasiparticles, we now proceed to derive the QBE 
that governs this distribution. We begin with the matrix 
form of Dyson's equation (Eq. [4]) In particular, we need 
to derive the equation of motion for G < . After a stan- 
dard short derivation, see for instance Refs. [l6| and [l3l . 
and working in the gradient expansion, we arrive at the 
expression, 



uj — Re E R , G < 

2m 



- [£<,ReG R ] 
S > G < -G > E < , (18) 



which describes the evolution of G < . Here [A, B] is a 
generalized Poisson bracket defined as 



[A,B] 



dA dB 
duj dt 



dA dB 
dt duj 



V r A ■ V k B 



V k A ■ V r B. 

(19) 



Note that in Eq. fTS! we have suppressed the variables 
(k, ui, r, t) for the self-energies and Green's functions. 

From standard perturbation theory working to one- 
loop order, the self-energies S < and S > are given by 



E 



r°° dv 


k X 


q 


JO T 


m 






k X 


q 


Jo * 


m 





ImD(q, v) [(n (v) + l)G<(k + q, uj + v) + n n {v)G< (k + q, uj - v)\ (20) 

InxD(q,z/) [n {v)G > {k + q, uj + v) + [n [v) + l)G>(k + q, uj - v)) , (21) 
I 



where for notation convenience, we have dropped the variables (r, t). Note that we have assumed that the 
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gauge bosons are always in local thermal equilibrium and 
that n (v) = l/(e f3v — 1), the standard boson equilibrium 
distribution function at temperature T. We therefore are 
studying the contributions to the transport coefficients 
arising from spinons. In Section VI we study the validity 
of this assumption for the particular case of the thermal 
conductivity. 

To derive the QBE for the generalized distribution 
function defined in Eq. [T6] we need to write Eq. [18] in 
terms of /(k,w,r,i). We integrate both sides of Eq. [TH] 
over the magnitude £ and rely on the assumption of the 
peaking as a function of £. From the Kramer's Kroenig 
relation, 



ReG fl (£,k,w) = -V 



dw' ImG fl K,k,w') 



7T tO — LO' 



(22) 



Since -2ImG i? = A, the condition that f #M = 1 i 
plies that 



|i R e G - = p/^^ = 0. (23) 

Z7T / Z7T tO — LO 



Therefore with the assumption that S < is independent 
of £, we can drop the second term on the LHS in Eq. [TH] 
and we are left with 



to — Re £ R , G < 

2m 



(24) 



We expand the remaining generalized Poisson bracket 
on the LHS and the QBE becomes, 



/*[( 



<9Re£ fl \ dG< dReT, R dG< 



dto J dt dt duo 

V r ReE fl • V k G< + V k [e k + Re£ H ) • VrG' 



(25) 



where the collision integral, I co ii, is defined below. Using 
our assumptions that the self-energies depend only on to 



even when the system is not in equilibrium and that G < 
remains a well-peaked function of £, we can perform the 
integration over £ and find that 



<9Re£ fl \ df dReZ R df 



V r ReS fl - V kF /+ 



dto ) dt dt dto 

V kF (e fe + Rc£ fl ) • V r / = Jcoii. 



(26) 



where / co n and £ now contain the generalized distri- 
bution function /(k, to,r,t) instead of the associated 
Green's functions. We have introduced the notation 
Vk F ff(k, to) which is defined as Vkff(k, to) evaluated at 
k = fcpk. Note that the term VkpReS^ • V r / is not zero 
in general because the self-energy can still depend on k, 
however, for this particular model the self-energy only 
depends on to so this term can be dropped. 

Eq. [221 is the full QBE for the generalized distribution 
function /(k, to, r, t). We see that despite the lack of a 
well-defined quasiparticle, this QBE looks very similar to 
the standard QBE derived for Landau Fermi liquid the- 
ory. Instead of the normal energy variable ek, the QBE 
now contains to which is independent of k. Moreover the 
QBE now involves renormalized time and energy deriva- 
tives of the distribution function^ 

We are interested in calculating the transport coeffi- 
cients in the linear response regime. In fact, we have al- 
ready assumed that the deviations from equilibrium are 
small, so that the generalized distribution functions are 
well-defined even out of equilibrium. We therefore pro- 
ceed to linearize the QBE of Eq. [26] By linearizing the 
QBE in a particular way, we are able to use a variational 
method to estimate the transport coefficients. Focusing 
on the RHS of Eq. (55] i- e. the collision integral, we in- 
troduce k' and to' as the energy and momentum of the 
intermediate spinon in the 1-loop self-energy diagram and 
find that the collision term, after the integration over £ 
from Eq. [24] becomes 



Icon = N{0) / dto'dk'dvdq ImD(q,z/) 



k'xq 



<5(fci?k' — kpk — q)x 

S(to' -u-u) [n (u)(l - f(k>, «'))/(£, to) - (no(i/) + 1)(1 - f(k u>))f(k', to') 

5(lo' -lo + v) [(noM + !)(!- f(k',to'))f(k,co) - n (v)f(k>, to')(l - /(k») 

I 



where N(0) = mfn is the density of states at the Fermi boson of energy v. Note that from Eq. [3j 
level for the up and down spins combined. Here we have 

rearranged the terms from E > and E < into the two pro- ^ V q 

cesses corresponding to absorbing and emitting a gauge InxD(q, v) = ^ 2j/ 2 + ^ 2 ' 



(27) 



(28) 
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where q — |q|. 

We define /(k,w) = f (uj) + 8f{k,Lu) 
some local equilibrium distribution, 

1 

/o(w,r, t) 



2 /3(r,t)(w-^(r,t)) + l' 



where /o(w) is 



(29) 



We re-iterate that all of the distribution functions both 
in and out of equilibrium are also functions of r and t. 
In particular the local equilibrium distribution fo(u>) can 
depend on space and time through the local temperature, 
(3(r,t). This local equilibrium solution has the property 
that it sets the collision terms on the RHS of the QBE 
to be zero/^ Thus from detailed balance, we derive the 
relations 



and 



no(l-/£)/o = (no + l)(l-/o)/o (30) 



(no + l)(l-/o)/o = no/5(l-/o), (31) 



for the 8{lu' — lu — v) and 8{lu' — lu + v) processes re- 
spectively. Note that for notational convenience we have 
introduced the definitions n = n a {u), f = fo(w) and 
/o = fo(u>'). We expect that for a fermion system the 
deviation from equilibrium 8f(k, lu) is sharply peaked 
around the Fermi surface, w = fj,. Therefore, we intro- 
duce a new function c/>(k, lu) defined by 



,dfa 
dou 

= /o + ^/o(l-/o). 



/ = fo - ■ 



(32) 
(33) 



Thus the function <fi is much smoother than the original 
8f. Writing the generalized distribution function in this 
way is also critical in deriving the variational method 
that we use to calculate the transport properties. 

Since we know that the local equilibrium distributions 
set the collision terms to be zero, we can expand out the 
distribution functions and using Eqs. 1301 and [3T1 find that 



Icon — 



J duj'd\t'dvd(\ ImZ?(q, v)x 



k'xq 



8(kpk' — kpk — q) x 
{S(lu' -lu- v)/3(<t> - 0> o /o(l - /o)+ 



(34) 



where = $>(k, lu) and <j)' 



(k'J 



III. FINITE TEMPERATURE 

Before proceeding to calculating the transport coeffi- 
cients using the linearized QBE, we must consider care- 
fully what happens at finite temperatures in this system. 
This problem was addressed by Kim et al. in Ref. H3 . 
We explain their argument here so that we can see how 
it effects our derivation of the QBE and later our deriva- 
tion of the transport coefficients. While in principle the 



derivations in section II are valid for both zero and fi- 
nite temperatures, in this system we need to take special 
care at finite temperatures because at any finite temper- 
ature T, the self-energy, ImSg (k, lu), is divergent even in 
equilibrium^ This is an infrared divergence and since it 
arises when the small q and lu limits have been treated ex- 
actly in cannot be alleviated through some sort of cutoff. 
One can show in fact that this divergence is an artifact 
of gauge invariance in the system. 

The divergence of the self-energy implies that the spec- 
tral weight A(k, lu) and the Green's functions are techni- 
cally not well defined at finite temperature, which seems 
to invalidate the basis of our derivation of the QBE. This 
problem can be addressed however by carefully consider- 
ing the source of this divergence. In particular, we see 
that the source of the divergence is the gauge field fluc- 
tuations where the energy carried by the gauge field v is 
such that v < T. We proceed then by breaking up the 
gauge field fluctuations into two pieces. Following Kim, 
we define a + (q, v) to be the fluctuations for v > T and 
a_ (q, v) to be the fluctuations for v < T. We then treat 
the a_ field as a vector potential which corresponds to a 
static applied random magnetic field b = V x a_ . Wc 
only consider the dynamics of the a + field. It is then 
understood that all equations need to be averaged over 
all possible configurations of this random magnetic field 
b . 

Kim et al. showed that after breaking up the gauge 
field one can regain the original QBE if the momen- 
tum k is shifted to k_ = k a and the self-energy 
is understood to only include fluctuations of the a + field. 
Since the self-energy now contains only fluctuations with 
v > T, it is no longer divergent. One then recovers the 
original QBE of Eq. [25] with an additional term on the 
LHS 



k_ 



• b_ x V k G < 



(35) 



corresponding to an applied random field b_. The origi- 
nal divergence of the full self-energy can be understood as 
a consequence of the non-gauge invariance of the original 
Green's function. We note that the extra term added to 
the QBE by breaking up the gauge fluctuations depends 
only on the gauge invariant quantity b_ and not on the 
potential a_ . 

We now proceed to consider the effect of this ap- 
plied random field on transport calculations using the 
linearized QBE that we have derived. In principle, we 
calculate the transport properties separately for each pos- 
sible configuration of the random field and then averages 
over them all. Because of this averaging, the field b_ 
cannot give rise to a linear response unlike an applied 
"electric field" or thermal gradient. It can however affect 
the transport coefficients of these quantities by contribut- 
ing an additional source of scattering. We can estimate 
the effect of this scattering by calculating a scattering 
rate due to a_ fluctuations. Kim et al. show that these 
fluctuations give rise to a scattering rate r _1 ~ T 4 / 3 . We 
can then use this estimate to calculate how this scatter- 
ing affects the low temperature forms of the transport 
coefficients. 
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IV. SPIN RESISTIVITY 



by 



The first transport property that we calculate is the 
spin resistivity. While less experimentally accessible than 
the thermal conductivity, it is technically simpler and 
thus serves to illustrate our method. We assume that 
there is some uniform applied force field F that couples 
linearly to the spinons. One possible way of realizing such 
a field is by applying a spatially varying magnetic field. 
Assuming that B is uniform along one direction and has 
a constant slope along the other direction, the spinons 
see a uniformly applied force field, F = VB. This field 
pushes up spins and down spins in opposite directions 
and thus leads to a net spin current, J s . 

This effective applied force field F couples to the 
spinons just as an electric field couples to an electron with 
unit charge. Thus we can borrow the result from a sys- 
tem of electrons in an applied electric field to determine 
the form of the driving term in the QBE. Following the 
work of Mahan in Ref. [lj, we find that in the presence 
of a scalar potential, the energy of the particle depends 
on its location, where as in steady state we expect the 
system to be spatially uniform. Mahan shows that this 
problem can be eliminated through a change of variables 
which has the secondary effect of changing the derivative 
V r to V r + F^j. This change of variables generates the 
driving term in the LHS of the QBE. 

We again linearize the QBE using the earlier definitions 
of fo(uj, r, t) and <fi(k, uj, r, t). Since the applied force field 
F is independent of position and time, we expect that 
that the steady state solutions for both fo(oj) and 0(k, uj) 
are independent of r and t. Thus returning to Eq. 1261 we 
see that the first two terms are zero. Using the spatial 
invariance of the distribution functions and the fact that 
the equilibrium distribution function is independent of k, 
we are left with 



F ■ kv F 



dfo 
duj 



coll- 



(36) 



This is the QBE that we consider for calculating the spin 
resistivity. Note that all of the terms containing the self- 
energy on the LHS have disappeared from the QBE. Thus 
the divergent effective mass does not even enter this cal- 
culation and should have no effect on the spin resistivity. 
We also note that Eq. [35] is almost identical to the lin- 
earized Boltzmann equation in an applied electric field 
with the change that k is now (k,ui). Thus we proceed 
by following Ref. [It] to derive a variational method for 
calculating the transport coefficients. 

We begin with defining the LHS of the linearized QBE, 
here Eq. [351 to be X(k, uj). We then rewrite the colli- 
sion term by defining a function P(k, uj, k', uj') that is the 
analogue of the equilibrium transition rate between the 
states (k, uj) and (k',uj'). Therefore the linearized QBE 
of Eq. [36] can be written as 

X(k, uj) = N(0) J dk'duj' (<j) - <j>') P(k, uj, k', uj'). (37) 

In order to proceed to derive the variational method we 
first define an inner product of functions of (k, uj) given 



(g,h) = N(0) / dkduj g{k, uj)h(k, uj) 



(38) 



We also define the operator V which takes the function 
(j) and returns the function given by 

T<t> = N(0) / dk'dw' {<t> - (j)') P{k, uj, k', uj'). (39) 



We note that the equilibrium transitional rate P is sym- 
metric in the exchange of (k, uj) with (k', uj'). It is easy to 
see that the operator V is linear and (</>, V<fi) > 0. From 
these properties, we can derive that the solution </> of any 
linearized QBE that can be written in the form of Eq. [37] 
is such that it minimizes the quantity^ 



A = 



(40) 



We now return to consider the specific case of the QBE 
given by Eq. 1361 Since the up and down spins contribute 
identically to the spin current J s and we are interested in 
ps only up to a numerical prefactor, we can consider the 
current due to the motion of only one type of spin which 
we denote by J. For a current density J the energy den- 
sity dissipation rate is p s 3 2 . This can in turn be related 
to the rate of entropy density production, so that we are 
left with 



TS 



PsJ 2 , 



(41) 



Looking at the definition of the spin current in terms 
of the generalized distribution function, we find that 



J = N(0) J dkdujv F kf(k,uj) 



(42) 



It is important to note that none of the renormalizations 
from the self-energy enter this expression. These renor- 
malization factors instead appear on the time derivative 
of the particle density. 12 This fact is critical for relating 
the variational principle to the transport coefficients and 
also helps explain why the diverging effective mass does 
make the transport coefficient ill-defined. 
Because of the k independence of /o(w), 



J = -iV(0) / dkduj v F k(j){k, uj) 



duj 



Thus if X is evaluated for unit applied force field, |F| 
we are left with 



,X) 



(43) 



(44) 



Following Ref. we now consider how the entropy 
is related to the distribution function /. In general the 
Boltzmann equation balances the rate of change of the 
distribution function due to diffusion, external fields and 
scattering. We can view the QBE of Eq. [26] in the same 
light. From statistical mechanics we know that the en- 
tropy of a system of fermions in equilibrium is given by 

S = -k B N(Q) j [/ In / + (!-/) ln(l - /)] dkduj. (45) 
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We assume that this formula holds for small deviations 
around equilibrium. Differentiating with respect to time 
and linearizing, the rate of entropy production is given 
by 



S = 



N(0) 
T 



(ftfdkdtu. 



(46) 



Thus from the collision term, we see that the rate of 
entropy production due to scattering is given by 



TS = (<t>,V<l>) 



(47) 



In steady state the macroscopic rate of entropy pro- 
duction, i.e. the Joule heating, is equal to the entropy 
production due to scattering. Thus we can solve Eq. [IT] 
for p s and using Eqs. |4"41 and Wf\ we find 



(48) 



Comparing Eq. @ 
that the solution 



i] to the earlier definition of A, we see 
that solves the linearized QBE of Eq. 
1361 minimizes p s . We thus proceed using the standard 
variational method. By making a reasonable ansatz for 
the deviation from local equilibrium <j>, we can calculate 
an estimate for the spin resistivity p s . 

Although the collision term includes the two distinct 
processes of absorbing and emitting a gauge boson, we 
note that we need only calculate p s for one of the two 
processes since they contribute identically to the spin re- 
sistivity. We consider the trial function 



tj>(k,u>) =77(k-F) 



(49) 



where F is a unit vector in the direction of the applied 
field F and rj, which has units of k F /to, is small. This par- 
ticular deviation from equilibrium (j> can be interpreted 
as a shift of the Fermi surface in the F direction which is 
a reasonable guess since the shifted Fermi surface is the 
ansatz used to derive the standard T 5 dependence of the 
low temperature resistivity in a metal. 

Both the numerator and the denominator of Eq. 05] 
contain factors of rf so they cancel and we drop this 
factor from the rest of our calculation. Plugging in the 
particular form of the trial function <f> from Eq. [49] into 
Eq. 221 it is straightforward to calculate the denominator 
which is independent of T and given by 



,X) 2 = (mv F ) 2 . 



(50) 



Looking back at the definition of the operator V given 
in Eq. [39] and using its symmetry properties, the nu- 
merator of Eq. 35] up to a numerical prcfactor is given 
by 

(0,7^) = J dLjdw'dvdqdkdk' \k' x q| 2 x 

InxD(q, V )tf - cj>') 2 f (oj)(l - /o(o/))x 

n M6(u}' -w-i/Wk' -k- -^-). (51) 

fcp 



We first perform the integration over ui and u/, using the 
result that 



dto du> S(u> — lu — i/)/o(cj)(1 — fo(u> )) 



p(i-e-*y 



(52) 



where z = 0v. We define 9 and 9' as the angles between 
k and u, and k' and u respectively. We then define a — 
9' — 0. The integration over q just enforces the condition 
from the momentum delta function, so that we are left 
with 

(4> ) v<f>) = p f dvdede' m2 f F - EiM - i m D( q (a),v)x 



1 — cos(a) 



(cos 8 — cos 6' 



0(1 



1 



(53) 



where we have defined the function q(a) which gives the 
magnitude of q for a particular angle a. We shift the 
integration over 9 and 6' to a = 9' — 8 and a' = 8' + 8. 
Integrating over a' gives 



V<f)) =0 I dvda m v F 



2 2 sin 2 (a) 



1 — cos(a) 



TmD(q(a), z/)x 



sin ( — )(27r — a + sun 



0(1 -e~ z )(e z - 1)' 
(54) 



We now assume that since ImD(g, v) is peaked for small 
q <IC /cf that we can take a to be small. Re-introducing 
q through a change of variables of a, we find 



b,Vcl>)=0 J dvdq^-ImD(q,v) — 



0(1 -e~ z )(e z - 1) 

(55) 

Inserting the RPA propagator from Eq. 1281 we arrive 



at 



,V4>) = 



yl/3 



xtkv0^ 




r 4/3 



(1 - e- z ) (e z - 1) 



y(q=2k F ) 



dy 



y 



1/3 



i + y 2 



(56) 



The lower limit of the integration over z is an artifact of 
our treatment of the divergence of the self-energy at finite 
temperatures as noted in section III. The self-energies in 
the QBE are taken to only include gauge fluctuations 
with v > fceT. We consider the effect of the low energy 
gauge fluctuations below. 

Looking at the integrand of the z integral in Eq. 1561 we 
see that it is sharply peaked for small z. Since y ~ 1/zwe 
can thus take the limit of the integration over y to be oo. 
The integrals over y and z are thus numerical prefactors 
of order unity that we ignore. Combining Eqs. 35] [50] 
and 1561 and reinserting the correct factors of h, we find 
that the spin resistivity is, up to a constant of order one, 



p s = h 



k B T 



4/3 



(57) 
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The final result for the spin resistivity takes exactly the 
form predicted by naive arguments assuming the exis- 
tence of quasiparticles and ignoring the effects of the mass 
renormalization. Thus the only effect of the gauge bosons 
is through the relaxation rate. However, we have now de- 
rived this result without assuming the existence of quasi- 
particles and we have shown that a possible divergence 
due to the effective mass does not enter the expression 
for the spin resistivity. 

Finally, we briefly discuss the effect of the low energy 
gauge fluctuations that cause the divergence of the self- 
energy at finite temperatures. As described in section 
III, these fluctuations enter the QBE as an applied static 
random magnetic field b_. Averaging over all possible 
configurations of this field does not lead to a linear re- 
sponse. Its only effect is through the transport scattering 
rate t_ ~ T~ 4//3 . Taking the result from the QBE deriva- 
tion, we assume that we can use the naive form for the 
spin resistivity ignoring the mass renormalization. Thus 
we expect that the low energy fluctuations in the self- 
energy loop do not change the overall scaling of p s and 
only enter through the numerical prefactor. 



V. THERMAL CONDUCTIVITY 

We now proceed to calculate the more experimentally 
relevant transport coefficient, the thermal conductivity. 
We consider the situation where a uniform thermal gra- 
dient is applied to the system giving rise to a heat cur- 
rent density U which is related to the thermal gradient 
through U = KV r T where k is the thermal conductivity. 
The heat current density U is defined in terms of the en- 
ergy current density Je and the particle current density 
J p , where 



U = J E - fiJ t 



(58) 



Note that the particle current is not equivalent to the 
spin current J s . Despite the fact that the quasiparticle 
is not well defined, these currents can be related to the 
generalized spinon distribution function f(k, w) through 

J E = N(0)^2 I ojv F kf a {uj,k)dkduj (59) 

J p = iV(0)^ j v F kf a (u),k)dkduj. (60) 

As in the case of the spin current, we again see that the 
renormalizations due to the self-energy do not enter the 
expressions for the currents and instead appear on the 
associated time derivatives of density or energy^ 
Thus the heat current density U is given by 



U = N(0) 



v-pk(u> — n)f a {uj, k)dkduj. 



(61) 



In the case of an applied thermal gradient, spin up and 
spin down react identically so we can drop the spin index 
for the remainder of the calculation. 



For a system in steady state in an applied thermal gra- 
dient, its clear that the local equilibrium fo(u>) defined in 
Eq. [29] depends on position r through a local tempera- 
ture /3(r) _1 . The distribution function however remains 
time independent. Thus the QBE for the system under 
these conditions is 



V r ReS K • V kF / + V kF e fc • V r / = /, 



coll- 



(62) 



We again linearize, assuming that / is expanded into 
/o and 4* as defined in Eq. [33] and that VT/T is small. 
Because fo(u>) is assumed to be k independent, the lin- 
earized QBE becomes 



Vk F e fe • V r /o = i, 



coll- 



(63) 



The earlier derivation of the variational principle re- 
mains true for Eq. [63] as well. Thus the quantity A, 
defined in Eq. [30] is still minimized by the solution 4> 
of the linearized QBE. Note that A is dependent on the 
particular QBE that we are considering. To proceed, we 
relate A for the linearized QBE of Eq. [63] to the ther- 
mal conductivity and then, in an identical way to the 
spin conductivity calculation, guess a reasonable <fr and 
calculate the associated thermal conductivity n. 

We again start with an expression for the rate of en- 
tropy density production^ For a given thermal current 
density, we have an entropy flux given by U/T. This 
leads to a rate of entropy density production 



s 



V(i) 



u 



u 2 

kT 2 ' 



(64) 



since U is constant throughout the sample. Note that we 
are assuming that there is no particle current due to the 
applied thermal gradient and thus no entropy production 
associated with a dissipation due to the presence of such a 
current. This assumption is valid for this system as long 
as there is no source of spinons. To see this we return to 
the Lagrangian of Eq. [1] Integrating out the gauge field 
a leads to the constraint that the spinon particle current 
density is conserved. Thus, the thermal resistivity W — 
is given by 



W = 



T 2 S 
U 2 ' 



(65) 



The derivation of S in terms of <p and V(f> from Eq. [47] 
is still valid here. Thus we only need to consider (cj), X) 
and see if it is related to the heat current density U. We 
again define the LHS of the linearized QBE, now Eq. [S3] 
to be X. Thus, 



, X) = AT(0) / dkduj(/)v F kV r fo 



(66) 



dfo 



= -N(0) / dkdu}(t>v F p(u-iJ,)-^-V T T-k. (67) 

J OLU 

Inserting the expression for / from Eq. [32] into the ex- 
pression for the heat current density U from Eq. I61[ 
gives an identical expression up to a factor of (3. We find 
that 



b,X) = -pV r T-lJ, 



(68) 
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and U 2 is thus equal to the denominator of A up to 
some temperature independent constant. Therefore the 
<j> that solves the linearized QBE of Eq. l63l minimizes the 
thermal resistivity W. Given an arbitrary <f>, W can be 
calculated, up to some numerical prefactor, with 



W 



T(</>,P<j>) 

W>,a} 2 ' 



(69) 



We need a trial function <fr that leads to a net flow of 
heat but no spinon particle current. A trial function that 
almost realizes this condition is 



(f> = r\ (to — /i)k ■ 



(70) 



where u is a unit vector in the direction of the heat cur- 
rent and r\ is small. This trial function can be interpreted 
as taking the original local temperature 0(r) and giving 
it a k dependence, 0(h) = + 77k • u. Physically at some 
point r, the spinons with momentum pointing in the di- 
rection that the temperature decreases are hotter than 
average and on the opposite side of the Fermi surface the 
spinons are colder than average. 

With the trial function (j) we now proceed to calcu- 
late the thermal resistivity W. First it is again clear 
that W has no rj dependence. Inserting the defini- 
tion of the trial function <j> from Eq. [70] we find that 
|U| ~ to?jf(&bT) 2 . In order to calculate (cf>,V4>}, we 
must recalculate (<j> — 4>') 2 for the new trial function. Us- 
ing the relations imposed by the energy and momentum 
delta functions, we find 

{<f> ~ <f>') 2 =v 2 (k • u) 2 + (w 7 - /i) 2 (q/^F • u) 2 + 

2z/(u;-/j)(k-u)(q/fc F -u). (71) 

Each of these terms can now be integrated individually 
in an analogous way to the calculation of the numerator 
of the spin resistivity. It is again necessary to use the fact 
that ImZ?(q, v) is peaked for small q in order to arrive at 
an analytic result. In evaluating these integrals, we need 
to calculate integrals of the form 

I{n)= f du>du>'5{u' -w-i/)/oH(l-/o(w'))(w-M) n , 



(72) 

for n = 0, 1,2. The integrals can be evaluated and we 
find 



1(0) = 



0(1 -e->)' 



1(2) 



20 2 (l-e- z Y 

TT 2 Z + Z 3 

30 3 (l-e~ z )' 



(73) 
(74) 
(75) 



Using these results, we can now calculate the integrals, 
(</>, V4>) for each of the three terms in Eq. [71] Dropping 
numerical prefactors, the first term in Eq. I71[ i.e. the 
v 2 term, contributes v\ (mfc?) 2 ^ 3 (fceT) 8 / 3 while the re- 
maining two terms both contribute 77i 4 / 3 fc p 2 ^ 3 '(k^T) 10 / 3 . 
Combining this result with our earlier result for (cf>, X) 2 



and plugging them into the definition of W from Eq. [69l 
we find that, restoring the correct factors of h, 



W 



k B e F 



1/3 



k B T 



1/3 



(76) 



Thus to leading order the low temperature thermal con- 
ductivity k per layer is 



a. - M ( jlS\ 

T~ h \k B Tj 



2/3 



(77) 



Thus we see that the QBE result for the thermal conduc- 
tivity agrees with the result we found by making the un- 
justified assumption that quasiparticles were well-defined 
and that we can ignore the effects of mass renormaliza- 
tion. Note that these are the same assumptions that gave 
the correct result for the spinon resistivity. 

With these assumptions we can show that k/T ~ te- 
We now use this result to consider the effects of the low 
energy a fluctuations. The applied static random field 
b_ again only enters the thermal conductivity through a 
relaxation rate since it cannot lead to a linear response. 
We therefore are interested in the contribution to the 
energy relaxation rate due to b_ fluctuations. This con- 
tribution goes as T~ 2 / 3 . Thus if we assume that the QBE 
result that k/T ~ te is valid , the low energy gauge fluc- 
tuations only change the numerical prefactor of k and 
not its temperature dependence. 

The form of the thermal conductivity from Eq. [77] 
is valid for clean systems only. In reality the divergent 
behavior of k/T at low T will be cutoff by the impurity 
scattering rate 1 /r . Again assuming that the QBE result 
is valid we thus find that the overall result for the thermal 
conductivity per layer in the presence of impurities is 



K 
T 



k 2 



k B T 



CF 



2/3 



mA 1 



(78) 



where A is the area of the layer. 



VI. THERMAL CONDUCTIVITY OF GAUGE 
BOSONS 

In all of the calculations up to this point we have as- 
sumed that the gauge bosons are always in local thermal 
equilibrium. Thus when we wrote down the QBE for 
the spinons, we assumed there was no deviation from 
the standard form for no(^). In this section we examine 
the validity of this approximation and in particular look 
at the corrections to the thermal conductivity that arise 
when this assumption is relaxed. 

As mentioned in the introduction, in the simplest 
derivation of the thermal conductivity, one finds that k 
is proportional to the heat capacity of the particles. It is 
known that for the model of Eqs. [1] and [3] the abundance 
of soft gauge field fluctuations produce a T 2 ! 3 contribu- 
tion to the specific heati^ As we showed in the previ- 
ous section the effective mass drops out of the expression 
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for the thermal conductivity due to the spinons. Ignor- 
ing the mass renormalization, the specific heat from the 
spinons is linear in T . Therefore, it is possible that at 
low temperatures, the largest contribution to the thermal 
conductivity comes from the gauge bosons. We now con- 
sider the effect of allowing them to deviate from thermal 
equilibrium. This effect is the equivalent of phonon-drag 
in the electron-phonon system. 

In order to describe the gauge bosons when they are 
not in thermal equilibrium, we derive the equation of 
motion for the gauge boson propagator that is analogous 
to Eq. [TH1 We then define a boson distribution function 
and derive an associated QBE. Thus we are left with two 
coupled differential equations that describe the system. 

Exactly as in section II, we begin by defining a num- 
ber of different gauge boson propagators, in analogy to 
the spinon propagators in Eqs. [6l through 1 1 1 1 except with 
t/j replaced by a. We denote these propagators with D 
rather than G. Note that this is consistent with the def- 
inition from Eq. [3] which we now associate with D . 
We also define an associated self-energy IT(q, v) which 
arises from the spinon bubble. In order to avoid dou- 
ble counting we consider the integration that lead to the 
gauge propagator in a renormalization group sense. In 
other words, we consider integrating out the high energy 
spinons to generate the gauge propagator and then con- 
sider coupling of the gauge propagator back to the low 
energy fermionic modes. 

Following the steps of section II, we consider the ex- 
pression 

D> (q, u)-D< (q, v) = -2 ImD R 

- 2lVq (79) 



2 „6 1 



7 ^ +XM 

It is clear that this expression is peaked at v — ^ q = 
(xd/t)? 3 - We can integrate v over the region of peaking 
and define the gauge boson distribution function as 



n(q., r, t) = J t^D< (q, u, r, t). 



(80) 



Note that unlike in the derivation of the generalized 
spinon distribution function, here we integrate over en- 
ergy just as one does in the standard Fermi liquid case. 

We now extended the assumption of the peaking of D 
to situations near thermal equilibrium as well. We see 
that despite the over-damped mode, the system looks 
as if it has quasiparticles at v = i/ q . This leaves us 
with the QBE for the gauge boson distribution function, 
analagous to Eq. [26] described by 



— + V q f„ • V r n = n< (n + 1) - U>n. 
at 
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We can now consider the two coupled QBEs of Eqs. 
[26l and [5T1 to calculate the thermal conductivity. Before 
proceeding, we linearize the gauge boson QBE. In anal- 
ogy to Eq. 1331 we define the deviation from equilibrium 
to be 



i(q) = n (u) - C(q) 



dnpiy) 
dv 



(82) 



noW+/3C(q)"oM(l + noM). (83) 



Note that in the above expressions v is technically f(q) 
since we have defined the gauge boson distribution func- 
tion as n(q, r, t). 

In principle, we would need to solve the coupled set of 
QBEs given by Eqs. [26] and [HTJ However, in the case 
of thermal conductivity there is an important simplifica- 
tion. We can calculate the thermal conductivity due to 
the spinons and gauge bosons independently, assuming 
in each case that the other type of excitations remain 
in thermal equilibrium. This is valid because in this as- 
sumption, the terms that are ignored can be shown to 
be smaller by a factor of k^T/e^ which is small for the 
temperatures under consideration in this system^ 

Thus to proceed we only need to calculate the addi- 
tional contribution to the thermal conductivity from the 
gauge bosons scattering off of spinons that can be as- 
sumed to be in thermal equilibrium. With this assump- 
tion one can evaluate II < and II < . Remembering that 
the gauge propagator Z?(q, v) of Eq. [3] was derived from 
the spinon bubble, we find that after linearizing, I co \\ for 
the boson QBE is given by 



I co i\ — J dcj'dk'dcjdk Im_D(q, w' — w)/?C(q)x 
noa-Zo^o^Fk'-^k-q). 



(84) 



We now make the ansatz £(q) = r/q ■ u. Exactly as 
before one can derive an expression for the thermal con- 
ductivity in terms of a variational method where 



1 



T 2 S 
U 2 ' 



(85) 



Just as in section V, the numerator can be related to 
the collision integral using Eq. 27] with 0(k, lu) replaced 
by £(q) and V now defined based on the gauge boson 
QBE. In order to use the variational method, we also 
must consider the left-hand side of the linearized form of 
the gauge boson QBE in Eq. [81] and relate it to U. 



(C,X) = J <iqC(q)v q • V r n K) 

= J dqC(q) (v q • VT) J»(q) 



(86) 
(87) 



VT 
~~T~ 



■ U 



since for the gauge bosons 



u = J dqn(q)^ q v q 



(89) 



Therefore we can use the variational method to calculate 
the gauge boson contribution to the thermal conductiv- 
ity. 

Returning to the denominator of Eq. [55] we see that 



<C(q),X) 2 ~T 2 



dn 
dT 



(90) 



Plugging in our variational ansatz for C(q)j one can show 
that 



T Z C 



(91) 
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where, as mentioned above the specific heat due to the 
gauge bosons C ~ T 2 / 3 . 

We now plug in the integral form of (£, V() into Eq. 
1551 giving 



1 



1 



Kg TC 2 



dujduj'dqdkdk'lmD(q,uj' - w)/?C 2 (q)x 
n Q (l- f )fo5(k F k' -k F k-q). (92) 



For the particular ansatz £ that we are considering, this 
integral is nearly identical to the calculation of the spin 
resistivity. After some work and borrowing the results 
from earlier calculations in this paper, we find 



1 

K a 



_!_T 4 / 3 ~ - 

TC 2 T' 



Therefore the total thermal conductivity is 



k ~ c s T 1/3 + c g T. 



(93) 



(94) 



At low temperatures the thermal conductivity is thus 
dominated by the contribution from the spinons and one 
can ignore the equivalent of the phonon-drag. 



VII. CONCLUSION 

Motivated by the organic compound k-(BEDT-TTF) 2 - 
Cu2(CN) 3 which is a quasi two-dimensional effectively 
isotropic spin 1/2 Heisenberg model on the triangular lat- 
tice, we studied a particular t/(l) spin liquid defined by 
Eqs. Q]and[3l that has been recently proposed to describe 
this system at low temperatures. From this model, a va- 
riety of properties that can be compared to experiment 
such as the specific heat and static spin susceptibility 
have been calculated. In this paper we have continued 
along these lines by calculating the transport properties 
of this spin liquid. 

In addition to the physical relevance of this calcula- 
tion, it is also interesting from a theoretical point of view 
because of problems that arise in deriving these coeffi- 
cients. In particular, there are three main issues that we 
had to circumvent. First, the renormalized spinons are 
not well-defined Landau quasiparticles because the one- 
loop correction to the spinon self-energy scales as u 2 / 3 . 
Second, the effective mass diverges at the Fermi surface. 
And finally, as an artifact of gauge-invariance, the self- 
energy is divergent at finite temperature. These issues 
clearly invalidate the results one finds from naive argu- 
ments using renormalized parameters. 

Despite these problems, we were able to proceed by 
deriving a QBE for this system. Starting from Dyson's 



equation, we showed that even though there is no well- 
defined concept of a quasiparticle, we can construct a 
generalized distribution function based on the peaking of 
the spectral weight as a function of instead of the 
standard u> peaking. Using this peaking, we derived the 
QBE for this system in terms of the generalized distribu- 
tion function. 

We then proceeded to linearize the spinon QBE and to 
calculate the transport properties. Applying the varia- 
tional method for solving the Boltzmann equation to the 
QBE, we solved for the transport coefficients. In partic- 
ular, we showed that the transport coefficients are well- 
defined despite the diverging effective mass. Moreover, 
we calculated the temperature dependence of the spin re- 
sistivity and thermal conductivity. Assuming the gauge 
bosons remain in thermal equilibrium, we found that at 
low temperatures, the spin resistivity goes as T 4 / 3 and 
the thermal conductivity goes as T 1 / 3 . In both cases, 
the result is identical to the one derived from naive ar- 
guments, where all the renormalizations due to the one 
loop corrections are ignored and the only effect of spinon- 
boson interaction is through a scattering rate 1/t. 

In the final section of this paper, we relaxed the as- 
sumption that the gauge bosons are in thermal equilib- 
rium. We showed that one can follow a similar method 
to the derivation of the spinon QBE and derive a QBE 
for the gauge bosons. In principle, we would then need to 
solve a coupled set of differential equations; however, for 
the specific case of the thermal conductivity we can de- 
couple the equations if we ignore terms of order k^T '/e-p. 
Under this assumption we showed that the gauge bosons 
contribution to the thermal conductivity is sub-dominant 
to the spinon contribution. 

Finally, we comment on the specific measurement 
of these transport coefficients in k-(BEDT-TTF) 2 - 
Cu2(CN)3. As mentioned in the beginning of section III, 
the spin resistivity is difficult to measure because of the 
problem of creating a source and drain of spinons. The 
thermal conductivity however is experimentally accessi- 
ble. Recent experimental and theoretical work on the 
specific heat and spin susceptibility have show evidence 
of spinon pairing at low temperatures (around 6K) in 
this system. Clearly the results in this work are invalid 
in the pairing regime; however, the exchange coupling is 
estimated to be J ~ 250K, so there is a wide region of 
temperatures where one could see evidence of the spin 
liquid state through the temperature dependence of the 
thermal conductivity. 
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